Adaptive Differential Feedback in Time- Varying 
Multiuser MIMO Channels 

Muhammad Nazmul Islam and Raviraj Adve* 
Dept. of Elec. and Comp. Eng., WINLAB, Rutgers University, Email: mnislam@winlab.rutgers.edu 
*Dept. of Elec. and Comp. Eng., University of Totonto, Email: rsadve@comm.utoronto.ca 



(n: 

O' 

(N; 
^— > . 

O: 
ov 



t/5 

o 



>: 

On 1 
O' 



X 



Abstract — In the context of a time-varying multiuser multiple- 
input-multiple-output (MIMO) system, we design recursive least 
squares based adaptive predictors and differential quantizers to 
minimize the sum mean squared error of the overall system. 
Using the fact that the scalar entries of the left singular matrix 
of a Gaussian MIMO channel becomes "almost" Gaussian 
distributed even for a small number of transmit antennas, 
we perform adaptive differential quantization of the relevant 
singular matrix entries. Compared to the algorithms in the 
existing differential feedback literature, our proposed quantizer 
provides three advantages: first, the controller parameters are 
flexible enough to adapt themselves to different vehicle speeds; 
second, the model is backward adaptive i.e., the base station and 
receiver can agree upon the predictor and variance estimator 
coefficients without explicit exchange of the parameters; third, 
it can accurately model the system even when the correlation 
between two successive channel samples becomes as low as 0.05. 
Our simulation results show that our proposed method can 
reduce the required feedback by several kilobits per second for 
vehicle speeds up to 20 km/h (channel tracker) and 10 km/h 
(singular vector tracker). The proposed system also outperforms 
a fixed quantizer, with same feedback overhead, in terms of bit 
error rate up to 30 km/h. 

Index Terms — Adaptive Differential Feedback, Scalar Quanti- 
zation, Multiuser MIMO Channels. 



I. Introduction 

The advantages of multiuser multiple-input-multiple-output 
(MIMO) has led to its inclusion in the standard proposals for 
fourth generation wireless systems, e.g., mobile WiMax (TJ. 
The best MIMO system performance can be achieved when 
channel state information (CSI) is available at the transmit- 
ter iff). In a frequency division duplexing system, downlink 
CSI needs to be estimated at the receiver, quantized and 
provided to the base station via an uplink feedback chan- 
nel. Recent work suggests that this might also be required 
in broadband time division duplex systems 0. Therefore, 
reducing feedback overhead while providing accurate CSI 
plays an important role in multiuser transceiver design. In the 
available literature, scalar quantization (SQ) |4j, 0, vector 
quantization (VQ) [6| and matrix quantization J7) have all 
been used to quantize CSI. SQ has been included in several 
standards (e.g. IEEE 802.1 In J8|) due to its linear compu- 
tational complexity. This paper focuses on the reduction of 
feedback overhead in time varying channels that employs SQ. 
We assume perfect channel estimation and delay-free noiseless 
feedback and focus on quantization only. 

The feedback overhead can be significantly reduced using 
differential feedback by exploiting temporal correlation of the 



channel ll^- lflTll . Most of these works model the channel 
as a first order Gauss-Markov process and use a fixed dif- 
ferential quantizer. The authors assume that the transmitter 
and receiver agree on the value of the parameters in the 
Markov chain. Since this assumption does not hold in non- 
stationary channels, there has been some research in adaptive 
delta-modulation (ADM) based feedback H, J5]. Both these 
works quantize the difference between the previous and current 
samples with a one-bit quantizer. However, the authors give 
suitable step size controller parameters for pedestrian veloci- 
ties (up to 4 km/h) only. 

The lack of flexibility of the available differential feed- 
back methods motivates us to investigate adaptive differen- 
tial feedback in time-varying multiuser channels. This paper 
makes the following two contributions: First, based on the 
linear least square (LLS) based adaptive differential speech 
quantizer model proposed by Stroh lfT2ll . we develop a 2-bit 
recursive least square (RLS) adaptive differential feedback in 
a time-varying environment. Second, we design RLS adaptive 
tracking of the singular vector entries of each users' channel 
matrix and show that, if the number of data streams is less 
than the total number of receive antennas, this method reduces 
feedback overhead. Both these methods can lead to reducing 
the required feedback overhead by several kBits/sec in modern 
wireless communication standards. 

Notation: Lower case, e.g., n or Lk, denote scalars while 
lower case bold face, e.g., h means a column vector. Upper 
case boldface, e.g., V denotes a matrix. The superscripts 
( ) T and (-) H denote the transpose and conjugate transpose 
operators respectively, tr [•] indicates the trace operator. E[-] 
denotes the expectation operator. I is reserved for the identity 
matrix. diag(x) denotes a diagonal matrix with non-zero 
entries taken from x. A(:, 1 : L) denotes the leftmost L 
columns of A. 

This paper is organized as follows: Section ITT1 develops the 
system model while Sections [III] and [IV] develop adaptive 
differential quantization for the channel entries and singular 
vectors respectively. Section [V] presents the simulation results. 
Finally, Section [VI] wraps up the paper. 

II. System Model 

We consider linear precoding with quantized channel knowl- 
edge in a multiuser MIMO environment where each users' 
vehicle is moving at an unknown speed. Specifically, we 
consider a single base station equipped with M transmit an- 
tennas and K independent users. User k has Nk antennas and 
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(a) Multiuser downlink 

Fig. 1. Block Diagram of Multiuser MIMO Downlink 

TABLE I 
(Channel Parameters) 



Parameter 


Value(units) 


Carrier Frequency (f c ) 


2.5 GHz 


Channel Sampling Rate (f s ) 


200 Hz 


Frame Duration (Tf r ) 


5 ms 
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Fig. 2. Block Diagram of the adaptive differential quantizer 



receives Lk data streams. Let L = J^k Lk and N — J^k Nk- 
Let U € C MxL and p £ 7Z L represent the beamformer 
matrix and power allocation vector respectively. P = diag(p). 
tr[P] < Pmax where P ma x is the total transmission power. 
The overall data vector is x = [x\,X2, ■ ■ ■ , xl\. The Nk x M 
block fading channel, Hf , between the base station (BS) and 
the user is assumed to be flat. The singular value decompo- 
sition of H fc is given by H k = A fc £ fc Bfc. S fe e K MxNk 
contains the singular values, global channel matrix is H H , 
with H= [Hi,...,H fc ]. 

In the downlink, user k receives 



yf L = H*UVPx + n fc , 



(1) 



where represents the additive white Gaussian noise at 
the receiver with E [nn fl ] = <r 2 lN k - To estimate its own 
transmitted symbols, from yf L , user k forms x& = V^y^ L . 
Let V be the N x L block diagonal global decoder matrix, 
V = diag (V 1 ,...,V K ). Overall 



x = V H H H U\/Px + V H n = F^UVPx 



V ff n 



(2) 



We define the MxL matrix F = HV with F = [Fi, . . . , F K ] . 
Ffc = HfeVfc. To ensure resolvability, L < M and Lk < Nk- 
In |[T3l . we showed that the sum mean squared error of the 
whole system can be written as, 



SMSE 



L 
i=l 



(3) 



- L - M + (a 2 + a%P max ) tr [J" 1 ] (4) 

J = FQF H + {a 2 + a e P m ax) 1m- o% denotes the quantiza- 
tion error variance of the feedback model and Q represents the 
virtual uplink power allocation matrix. d4]> is a nonincreasing 
function of a 2 E . Our objective is to minimize <j 2 E by designing 
a backward adaptive differential feedback system. 

Channel Mobility: The channels are temporally corre- 
lated and assumed to follow a modified version of Jakes' 
model lfl4l . Here, channel parameters are selected (as in 
TablelJi to represent typical values for the WiMax standard (TJ. 

Feedback Model: We use two different feedback methods. 
In the first method, the full channel matrix is quantized and 



sent back to the base station (BS). The receivers expend 2 bits 
on the differential quantization of the real and imaginary parts 
of each scalar channel entry based on minimum Euclidean 
distance. The BS uses the following channel model: 



H = H + H 



(5) 



Here H and H denote the quantized channel and error in 
channel feedback respectively. This channel model is used 
to find the optimal F through iteration [15|. Also, a 2 E = 



E 



\h-hf 



h and h denote the original and quantized 
channel entry. 

In the second method, the receivers use, as Vfc, the Lk 
right singular vectors corresponding to the maximum singular 
values of Hfe. So, Vfc = Bfc(:, 1 : Lk) Therefore, F^ = 
Afc(:,l : Lk) x Sfc(:, 1 : Lk). The receivers perform 2 bits 
adaptive differential quantization of each real and imaginary 
scalar entry of Afc(:, 1 : Lk) and 2 bit fixed quantization of 
the entries of Sfc(:, 1 : Lk). The BS assumes the following 
model, 

F = F + F (6) 
Here, F and F represent the quantized effective channel and 
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error in the feedback respectively. Here, a 2 E 

f and / denote the scalar entries of F and F respectively. 

The linear precoding algorithms of these two feedback 
models can be found in |15| and |fl6l respectively. We do 
not include those here for brevity. 

III. Quantization of Channel Entries 

Stroh [12] proposed the differential quantizer model, shown 
in Fig. |2] for speech quantization. We use it to perform 
adaptive differential channel quantization. The left and right 
sides of the channel block are located at the receiver and base 
station respectively. A unit variance Gaussian quantizer ifTTl 
is used in the quantizer block. Let h n and h n represent the 
original and quantized channel parameter at the n th instant. 
d n is the difference between h n and predicted channel h n . g n 
normalizes the variance of the difference signal, i.e., avoids 
granular noise and overloading. Thus, g n enables the quantizer 
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TABLE II 
(Design Parameters) 



TABLE III 

CODEBOOK OF SCALAR ENTRIES OF LEFT SINGULAR MATRIX 



Parameter 


A 


k 2 


kj 


T 


L P 


Value 


0.98 


0.9 


1.1 


2 


100 



block to adapt to different speeds. d n = d n + qn n where qn n 
is the quantization noise. 

Using the symmetry of the receiver and BS, h n — h n + 
qn n [12], [ 1 8 1 . Here, the adaptor block controls the predictor 
coefficients and the variance estimator block estimates g n . 
The predictor coefficients and variance estimator parameters 
depend on h„ and d„, rather than on h„ and d„. Therefore, 
unlike the differential feedback model proposed in l9l- lfTT1 . 
the BS can reproduce the predictor and variance estimator pa- 
rameters without the explicit transmission of the coefficients. 
Since the real and imaginary part are quantized separately, 
g n = 2a%. 



A. Parameter Selection of the adaptive quantizer 

1 ) LLS Based Predictor: Stroh proposed the following LLS 
based predictor in the speech quantizer |fl~2], 



T 

w j,nh n -j 



(7) 



Here, T is the predictor order and Wj jTl is the j weight 
coefficient at the n th time instant. The predictor coefficients 
are computed to minimize the mean squared error 
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(8) 



Here, L p is the learning period. e is the fitting or average 
prediction error. The weights are calculated through Weiner 
filtering |fl9l . 

2) RLS Based Predictor: LLS predictors perform close 
to ideal Weiner filter predictors in terms of quantization 
error reduction. However, reducing the steady state error to 
acceptable levels requires increasing the learning period and 
an attendant increase of transient time. This motivates us to 
design a recursive least square based backward predictor 1 19 1. 
The weights w(n) are calculated by solving the Weiner-Hopf 
equation, $(n)w(n) = ip(n); w(n) is the T x 1 vector of 
predictor coefficients, at the n time instant. <fr(n) e lZ TxT 
and ip(n) £ lZ Txl are given by, 

*(n) = A*(n 1) + h(n- l)h H (n ~ 1) (9) 
tp(n) = \>if)(n- 1) +h(n- l)d H (n- 1) (10) 

Here, h(i) = hi, - ■ ■ , h^T+i and A is the memory factor. 



M = 2 


M = 3 


M = 4 


M = 8 


Standard Gaussian 


-1.34 


-1.40 


-1.43 


-1.48 


-1.51 


-0.43 


-0.44 


-0.45 


-0.45 


0.45 


0.43 


0.44 


0.45 


0.45 


0.45 


1.34 


1.40 


1.43 


1.48 


1.51 



3) RLS Variance Estimators: For the RLS variance estima- 
tor, g n is calculated as [19|, 



ki J (1 - fc 2 ) I k 2 v n -i + d 



(ID 



(12) 



Here, k\ and k 2 denote the bias compensator and memory 
factor respectively in RLS variance estimator. Table [TT] lists 
the parameters of the adaptive differential quantizer. The 
parameters were chosen via numerical simulations. 

IV. Quantization of singular vectors 

The scalar entries of Afc(:, 1 : Lk), the left singular vector 
matrix of H, can be adaptively differentially quantized using 
the same model shown in Fig. [2] Note that both the adaptive 
predictors proposed in the previous section do not assume 
any particular model of the signal; they try to find the "best" 
predicted value based on the past observations. However, the 
quantizer in the proposed adaptive differential feedback model 
assumes a Gaussian distributed input. Therefore, if we can 
show the entries of to be approximately Gaussian, the 
model of Fig. [2] can be readily applied to track Afc |20|. 

The matrix of singular vectors of a rectangular Gaussian 
matrix is called a Haar matrix l20l . 

Lemma 1: If A k is a C MxM Haar matrix, E [lA^I 2 ] = -jj 
, 1 < i,j < M. 

Proof: See EQl. _ 

Lemma 2: The probability distribution of \f~M times the 
Haar matrix A^, approaches the standard complex Gaussian 
measure as M — > oo. 

Proof: See 4.2.11 of GUI 

In practice, the entries of the Haar matrix approach a 
Gaussian random variable for small values of M. To show 
this, we set Nk = 2 and choose different numbers of transmit 
antennas, M. We generate 10 5 random Gaussian distributed 
channels, H/j e C MxNh , and find the left singular matrix 
of Afe e C Mx . We randomly pick different entries of 
Afe. After normalizing the samples using Lemma 1, we find 
the 2 bit codebook of the collected samples using k-means 
clustering (21] . In Table [Til] we compare the codebook with 
that of a unit variance 2-bit standard Gaussian quantizer ifTTIl . 
In Table [TTTJ "Af = 2" stands for the 4 level normalized 
codebook, based on the scalar entries of the left singular 
matrix with 2 transmit antennas. The table shows that even 
for small number of transmit antennas (e.g., 3), the probability 
distribution of the normalized scalar entries of the Haar matrix 
resembles the Gaussian distribution. 
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The degrees of freedom of the Haar unitary matrix is less 
than the total number of real and imaginary entries. The 
minimum number of parameters to represent the Haar matrix 
can be extracted through Givens' rotations [4|. An adaptive 
controller for tracking Givens' rotated parameters have only 
been provided for pedestrian velocities j4|. The phases and 
Givens' rotated angles are not Gaussian distributed and least 
squares based predictors are not optimum to track these param- 
eters. Therefore, we stick to our proposed adaptive differential 
quantization policy. This ensures greater adaptability of our 
model at the cost of slightly higher feedback overhead. 

A. Adaptive Differential Quantization of singular value 

The distribution of the square of the singular values of a 
Gaussian channel i.e., the eigenvalues of the Wishart matrix 
can be found in 11221 . We perform fixed 2 bit quantization 
of the singular values using this distribution and the standard 
Lloyd-Max quantizer. 

In this approach, receivers feed back F = HV to the 
BS, instead of providing H. The dimensionality of H and 
F are M x N and M x L respectively. Thus, singular vector 
quantization saves feedback overhead as long as L < N. 

B. Discussion 

Sections Hill and ITVl show that, since we assumed the channel 
to be Gaussian and the difference of two correlated Gaussian 
random variables leads to another Gaussian random variable, 
the model shown in Fig. [2] provides great flexibility and can 
hold for different vehicle speeds. To the best of our knowledge, 
this is the only work in adaptive differential limited feedback 
literature, which can provide both the following advantages: 

1) Unlike the Gauss-Markov models of ll9l- lfTTl . our model 
works when the normalized autocorrelation between 
successive channel samples drops below 0.5. 

2) Unlike the feedback model proposed by 0J, 0, the 
controlling parameters of the predictor and variance es- 
timator in our model do not depend on the knowledge of 
the correlation between two successive channel samples. 

V. Numerical Results 

The channels were generated using the channel model 
of lfT4ll for each speed. Figure [3] shows the quantization 



error variance (er^ = E 



\h — h\ 2 ) of different proposed 
methods. The quantization error variance plot of a 2-bit and 
3-bit unit variance fixed Gaussian quantizer were plotted using 
the standard values (0.1175 and 0.0345 respectively) [17]. 
The figure shows that the RLS and LLS predictors perform 
very similarly in terms of quantization error reduction. The 
performance of the ideal differential quantizer and predictor 
is governed by 

<4 = (13) 
a% ln = 0.1175^ (14) 

Here, ( fT~3b follows the minimum error surface of an ideal 
Weiner filter [19] and ( fT4b follows the quantization error 
associated with a 2-bit, unit-variance, Gaussian quantizer ifTTl . 
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Fig. 3. Comparison of of differential feedback with fixed feedback 
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Fig. 4. Comparison of the transient time of RLS and LLS based feedback 
at 21.6 km/hr 



The RLS adaptive differential quantizer's performance be- 
comes inferior as the vehicle velocity exceeds 32 km/h. 
Using the parameters from Table |I] this speed corresponds 
to a maximum normalized correlation of 0.0255 between two 
successive channel samples. Therefore, our proposed adaptive 
differential feedback outperforms fixed quantization as long as 
the normalized channel autocorrelation remains positive. 

Fig. H] shows that the transient time of RLS adaptive differ- 
ential feedback is much smaller than the LLS one. The simu- 
lation was performed at 21.6 km/hr. The average quantization 
error was calculated at every iteration. The error variance of 
the proposed quantizer converges close to its final value within 
20 iterations, i.e., 100 ms. Thus, the proposed differential 
quantizer can adapt itself in real time with reasonable changes 
in the mobile velocity. Previous works in differential feedback 
literature have either focused on stationary channels with 
fixed mobile velocity [9|-[H] or non-stationary channels with 
pedestrian velocity |4), 0. Our proposed model are suitable 
for vehicles whose velocity can increase up to 30 km/hr. 

Figure[5]shows the average bit error rate (BER) performance 
of different feedback models at different speeds. We also 
represent the respective feedback overheads in terms of kilo 
bit per second (kB/s). We used the system model of Section 
II, quadrature phase shift keying modulation and the linear 
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Fig. 5. BER and overhead of different feedback methods, M = 4, K = 
2, Ni = 4, 7V 2 = 4, Lx = 2, L 2 = 2 

transceiver design algorithms of |15j and lfl6ll to simulate 
the performance of channel and singular-matrix quantization. 
Here, "adap chan" and "adap eig" denote channel quantization 
and singular vector quantization respectively. Figure [5] shows 
that the 2-bit adaptive channel entry feedback outperforms 3- 
bit fixed feedback and performs very close to the full channel 
knowledge feedback scenario at 1 1 km/h . Even at a high speed 
of 30 km/h (corresponds to a normalized autocorrelation of 
0.1 with a degree arrival angle [23]), the proposed adaptive 
feedback reduces the BER by a factor of 2, with respect to 
2-bit fixed feedback per channel entry. 

In Fig. |5j "2 bit adap eig" indicates use of 2 bits to quantize 
each of the real and imaginary parts of in an adaptive 
differential manner. Since, we assumed N = 8 and L = 4 
in our simulation, spending 2 bits per scalar singular matrix 
entry is equivalent to spending 1 bit per real and imaginary 
scalar component of the channel. At low speeds like 1 1 km/h, 
the singular matrix entry quantizer performs approximately 
as well as the 2-bit fixed quantizer and reduces the feedback 
overhead by a factor of 2 for almost same BER. Thus both the 
adaptive differential feedback methods save 1 bit per real and 
imaginary entry of the channel matrix at low speed (20 km/h 
for the channel tracker and 8-9 km/h for the singular matrix 
tracker). This leads to a saving of 2MNF S bits in feedback 
overhead per second. Using Table U parameters, the proposed 
systems provide a feedback reduction of 12.8 kBit/sec. 

VI. Conclusions 

In this paper, we developed adaptive differential scalar quan- 
tization based limited feedback in a time varying multiuser 
MIMO channel. The key contribution is the development of a 
differential feedback system that tracks the channel variations 
without a priori knowledge of the correlation across time, 
especially the speed of the vehicle. 

We proposed two methods of adaptive differential feedback. 
First, we developed 2-bit adaptive differential quantization of 
each scalar real and imaginary entry of the MIMO channel. 
Second, we developed 2-bit adaptive differential quantization 



of each scalar real and imaginary entry of the singular vectors. 
Both these methods were shown to significantly reduce the 
feedback overhead - by as much as 12.8kBit/sec. Our proposed 
adaptive differential quantizer model was shown to be flexible 
enough to adapt to different speed of the vehicles. 

It is worth emphasizing that one issue not addressed here is 
adaptive tracking of the channel gains. When different vehicles 
are located at different distances from the base station, we 
believe that an adaptive tracking of the channel gain would 
outperform a fixed channel gain quantizer. 
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